Cross-dimensional relaxation in Bose-Fermi mixtures 
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We consider the equilibration rate for fermions in Bose-Fermi mixtures undergoing cross- 
dimensional rethermalization. Classical Monte Carlo simulations of the relaxation process are per- 
formed over a wide range of parameters, focusing on the effects of the mass difference between 
species and the degree of initial departure from equilibrium. A simple analysis based on Enskog's 
equation is developed and shown to be accurate over a variety of different parameter regimes. This 
allows predictions for mixtures of commonly used alkali atoms. 
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I. INTRODUCTION 

In the study of ultracold dilute atomic gases, the in- 
terparticle interactions are characterized by the s-wave 
scattering length a. Because of the low temperatures 
achieved, and the lack of long range or anisotropic in- 
teractions, the value of a determines a wide variety of 
equilibrium and dynamical properties of quantum degen- 
erate gases. Furthermore the parameter a can be tuned 
in many systems, from — oo to oo by means of Feshbach 
resonances pj, allowing the experimenter access to any 
desired interaction strength. This potential for real-time 
control over interactions is a unique feature of ultracold 
gas experiments. Additionally the efficiency of evapora- 
tive cooling relies on a large elastic collision rate, which 
is proportional to a 2 . For these reasons, the ability to 
accurately determine the scattering properties of dilute 
ultracold gases becomes essential for quantum gas exper- 
iments. 

In their pioneering work with ultracold 133 Cs atoms, 
Monroe and co-workers showed that a relatively sim- 
ple rethermalization measurement starting with a non- 
equilibrium gas could provide a determination of the elas- 
tic collision cross-section 2], which is equal to 87ra 2 for 
identical non-condensed bosons. The rethermalization 
rate for a gas in the so-called collisionless regime (de- 
fined by a collision rate per particle much lower than 
the harmonic frequencies of the trapping potential), was 
measured by selectively removing energy from the gas in 
one spatial dimension and watching the subsequent cross- 
dimensional rethermalization. The relaxation in such ex- 
periments is driven by elastic collisions, and an analysis 
based on Enskog's equation shows that the rate of relax- 
ation is proportional to the mean rate of collisions 0, Q , 



Trclax = - (n> a (v rcl ] 
a 



(1) 



where n is the number density of the gas, a is the elas- 
tic collision cross-section, and v re \ is the relative collision 
speed; brackets (•) denote a thermal average, and we have 
used the assumptions of energy- independent s-wave colli- 
sions and Boltzmann statistics to write the mean collision 
rate (r co n) = (n<Jv IC \) — (n}a(v re \). The constant of 
proportionality a, defined as the ratio of collision and 
relaxation rates, reflects the mean number of collisions 
per particle required for rethermalization. A variety of 
numerical and analytical studies have found that a is 
between 2.5 and 2.7 [IHUH^. 

We have recently extended the method of cross- 
dimensional relaxation to probe the s-wave scattering 
length asF between bosonic and fermionic atoms in a 
Bose-Fermi mixture 8]. In that work, we compared 
single-species boson rethermalization to relaxation of 
fermions in the presence of bosons in order to elim- 
inate the systematic calibration error associated with 
atom number that typically dominates cross-dimensional 
rethermalization measurements. Since collisions between 
spin-polarized ultracold fermions are forbidden by the 
Pauli exclusion principle Q, their rethermalization pro- 
ceeds only through collisions with the bosons. The col- 
lision rate per fermion in the mixture therefore depends 
only on the number of bosons, allowing us to write the 
mean relaxation rate per fermion in analogy with the 
single-species case, 



I> = \ (n B ) f bf (vbf) 
P 



(2) 
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Here (tib) is the equilibrium density of the bosons, av- 
eraged over the fermion distribution, gbf is the boson- 
fermion elastic cross section, and (vbf) — (8 ksT j-n\x) x l' 2 
is the thermally averaged relative collision speed 
between bosons and fermions with reduced mass 
fx = mp tub I (nip + ms) and temperature T {k B is 
Boltzmann's constant). The cross-section is related to 
the interspecies scattering length by ob f — 4iva BF . The 
constant of proportionality [3 reflects the mean num- 
ber of collisions per fermion needed for rethermalization. 
Knowledge of a and (3 was essential in achieving both the 
accuracy and precision of the measurement in Ref . • 
In this work, we study the general dependence of (3 on 
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the masses of the fermions and bosons by means of de- 
tailed classical Monte Carlo simulations of the relaxation 
process. We additionally address the effect due to the 
finite initial departure from equilibrium. Understand- 
ing this effect is a necessary component of the analysis 
of both the simulations and the measurements, where 
a larger initial perturbation improves the signal-to-noise 
ratio. We further develop a classical kinetic model based 
on Enskog's equation for the case of Bose- Fermi mixtures 
that reproduces the behavior of the simulations. The 
analysis and simulations show the following: (i) Equa- 
tion J2J is valid for a wide range of parameters relevant 
to current experiments, (ii) the difference in mass be- 
tween the bosons and fermions can lead to a ~ 5 times 
difference in (3 between light and heavy fermions in mix- 
tures of experimental interest, and (hi) the size of the 
initial perturbation must be taken into account to fully 
understand the results. 

The remainder of the paper is organized as follows. In 
Sec. |n]we describe the details of the Monte Carlo sim- 
ulations and use the simulations to verify the validity of 
Eq.(|2J). In Sec. IHII the mass dependence of the relaxation 
is investigated using the simulations and compared to a 
classical kinetic theory that is developed in some detail. 
It is shown that the model reproduces the behavior of 
the simulations over a wide range of masses. In Sec. IIVI 
the importance of the size of the initial perturbation is 
addressed. Finally, we conclude in Sec. [V] with a dis- 
cussion of possible extensions of the model and further 
applications of the Monte Carlo simulations. 



II. MONTE CARLO SIMULATION 

Implicit in the definition of (3 in Eq. J3J) are the assump- 
tions that the energy anisotropy undergoes exponential 
decay ,5] and that (3 depends only on intrinsic proper- 
ties of the particles under study, such as the mass and 
quantum statistics, and not the bulk properties of the 
gases, such as the temperature or atom numbers (assum- 
ing one remains always in the collisionless regime) . These 
assumptions were verified by running Monte Carlo simu- 
lations of the relaxation. 

In the experiment considered here, we start with a 
Bose-Fermi mixture that is in thermal equilibrium in a 
cylindrically symmetric harmonic trap and increase the 
radial trapping frequency by a factor of fl to produce 
the initial energy anisotropy in the system; the axial fre- 
quency is unchanged. The change in trap strength is per- 
formed slowly with respect to the radial trap periods but 
quickly compared to the thermal relaxation time. The 
increase in radial energies E x ^ y is then the same for each 
species and is equal to the fractional increase in the trap- 
ping frequencies. Therefore the ratio of radial and axial 
energies immediately before rethermalization is simply 



E 



x.y 



E z 



= n 



t=0 



Note that for gases prepared in this manner, E x « E y 
at all times. Since the gases rethermalize together, 
they reach the same final temperature. In the classical 
gas limit, the final temperature is given by (Too/To) = 
(1 + 2il)/3, where To is the temperature before com- 
pression. It is important to note that for this type of 
rethermalization, there is no net transfer of energy be- 
tween species during the relaxation process; this fact is 
essential for the validity of Eq. J5J ■ 

An outline of the Monte Carlo simulation is as follows. 
An ensemble consisting of Nb bosons and Np fermions is 
prepared in a cylindrically symmetric harmonic trapping 
potential. The ensemble is initialized by assigning a ran- 
dom position and velocity vector to each particle from 
separable Gaussian distributions. The distributions are 
scaled such that there is a factor of Vl imbalance between 
the mean energies per particle (both kinetic and poten- 
tial) in the transverse {x,y) and axial (z) dimensions. 
The initial energy in a given dimension, however, is the 
same for both species. The positions and velocities are 
evolved for some small time step At according to New- 
ton's laws in the trap, and then collisions are considered. 

Only boson-boson and boson-fermion collisions are al- 
lowed, in accordance with the Pauli exclusion principle 
for a system with spin-polarized fermions. If two particles 
are found within a critical distance r c of each other, the 
pair is given a chance to collide. The collision probability 
is given by P co \\ = &m v ic i At/V c , where ati is the collision 
cross section between particles k and I, v Te \ = \i)k — vi\ is 
the relative collision speed, and V c — Airr^/S is the vol- 
ume of the sphere containing the colliding atoms. If P co ii 
is greater than a uniformly distributed random number 
between zero and one, then the relative velocity vector 
is rotated into a random direction (conserving the total 
momentum and energy) to represent an s-wave collision. 
Finally, after all possible collision pairs have been con- 
sidered, the mean energies of each species in each of the 
three cartesian directions are recorded, and a new time 
step proceeds. 

In agreement with Ref.0 we find it essential for ob- 
taining robust results to set At and r c such that the 
ensemble average collision probability (P C oii) and occu- 
pation of the fictional spheres surrounding the atoms 
are both well below 10%. Typical values for our sim- 
ulations are At ~ 1CU 3 to 10~ 2 times the collision time 
(ncrw ro i) _1 , and r c ~ 1CU 3 to 1CU 2 times the root-mean- 
squared (rms) cloud radius. Since the motion between 
collisions is known analytically, it is not necessary to keep 
At small compared to the trap period. Finally, we note 
that we have not observed any double-counting of colli- 
sions in consecutive time steps. This is because our time 
step At typically satisfies w > r c , where v is the char- 
acteristic speed of the particles. 

A typical relaxation curve obtained from the simula- 
tion is shown in Fig. ^ The relaxation rate of the 
energy anisotropy of the fermions is determined by fit- 
ting the ratio of radial and axial energies to the ratio of 
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FIG. 1: Typical Monte Carlo simulation of cross-dimensional 
relaxation of fermions in the presence of bosons. The mean 
energy per fermion Ei in each direction is shown versus time 
during the equilibration process. The curves represent 5 x 10 4 
fermions in the presence of 5 x 10 4 bosons. The calculation 
used 40 K ( 87 Rb) as the fermion (boson), and the trapping 
frequencies, temperature, and initial perturbation were taken 
from the experiment in Ref. ||. A cross-species scattering 
length of \cibf\ = 235 ao, with ao the Bohr radius, was as- 
sumed. The calculated mean time between collisions is 0.23 s. 

decaying exponentials, 

1 E x + E y _ 1 + f exp(-I> t) 

2 E z ~ 1 - 2eexp{-T F t) 

Here, Ei is the mean total energy per particle of the 
fermions in the i th direction, and e = (f2 — 1)/(1 + 20) 
characterizes the amount of initial perturbation from 
equilibrium. 

The validity of Eq.Q was tested by running Monte 
Carlo simulations with varying temperatures, numbers 
of atoms, trapping frequencies, and interspecies cross- 
sections ubf- As an example, the dependence of Tp 
on the number of bosons (Nb) and fermions (Np) in 
the mixture is shown in Fig. [21 The results show that 
T^ varies linearly with Nb, while it is constant when 
changing Np over the same range. Although these sim- 
ulations assumed \(1bf\ = 235ao and clbb = 98.98ao 
[To|. results similar to those in Fig. [2] were obtained for 
o-bf = 50 ao, corresponding to a ratio of cross-sections 
(&bb/vbf) = 0.35 and 7.8, respectively. 

III. MASS DEPENDENCE — CLASSICAL 
KINETIC MODEL 

The Monte Carlo simulation is a powerful tool for an- 
alyzing the behavior of relaxing Bose-Fermi mixtures, 
but it requires a considerable amount of computing re- 
sources. For example, one of our simulations using 10 4 



FIG. 2: Fermion rethermalization rate Ff as a function of 
Nb and Nf- These calculations used the parameters from 
the 87 Rb- 40 K experiment in Ref.Q, with \clbf\ = 235 ao. 



fermions and 10 4 bosons with 10 4 time steps took over 
8 hours on a Unix system with a 2.4 GHz processor and 
4 GB of RAM. The computation time, which is domi- 
nated by the pair wise search for collision partners, scales 
roughly as the product of the number of time steps and 
Nb {Nf + Nb/2), so that simulations become impracti- 
cal for N > 10 5 . Furthermore, it is desirable to build up 
a physical picture of the relaxation that is not provided 
by the simulations. For these reasons, we now consider 
an analytic model of the rethermalization. 

One new degree of freedom for two-species rethermal- 
ization is the appearance of a second mass. In order to 
see how the number of collisions per fermion needed for 
rethermalization depends on the masses of the particles 
involved, we first consider the two limiting cases of very 
light and very heavy fermions. Recall that energy is re- 
distributed by means of s-wave collisions, which random- 
ize the relative velocities of the colliding particles. For 
very light fermions, a single collision should be nearly 
enough to redistribute the energy. If the fermions are far 
heavier than the bosons, however, it should take many 
collision times to redistribute the energy, since a single 
collision has little effect on the motion of the heavy parti- 
cle. We therefore expect the number of collisions needed 
for equilibration to be a monotonically increasing func- 
tion of the normalized fermion mass i] = mp/(mp + rriB)- 
Note that this is in contrast to the behavior of two gases 
initially at different (isotropic) temperatures and subse- 
quently brought into contact for rethermalization. In the 
latter case, one expects /3 _1 oc 4 77 (1 — rj) [llj, which is 
most efficient for 77 = 1/2. The results of our Monte Carlo 
simulations is shown as the solid points in Fig. [3] As ex- 
pected, (3 increases smoothly with the fermion mass. The 
open points and solid lines are predictions from a simple 
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1/2 

i/ = m_/(m_ + m.) 



FIG. 3: Dependence of the mean number of collisions 
per fermion needed for rethermalization on the normalized 
fermion mass. Solid points (•) are the results of the Monte 
Carlo simulations, and open points (o) are from exponential 
fits to our model from Eq.©. The solid lines are the bounding 
predictions from the classical kinetic model. The number of 
atoms is Nb = Nf = 10 4 , the initial perturbation is Q = 1.64, 
and the ratio of cross sections is (obb /<tbf) = 0.35. 



kinetic model, which we now discuss in detail. 

Our analysis of the rethermalization process is based 
on Enskog's equation, which is equivalent to the Boltz- 
mann transport equation Q. This treatment gives the 
rate of change of the ensemble average of any function of 
the boson and fermion positions and velocities, usually 
denoted x(%b, v*b] %f, &f)> by 



_X) = cr{n u re i 



(3) 



where we have again used the assumption of energy- 
independent s-wave scattering to separate the collision 
cross-section a from the ensemble average. The quantity 
A\ is just the change in x due to a single collision. Here 
we consider 



Xi — E\ x — Ex* 



(4) 



where we are focusing for now on generic particles of 
type 1 colliding with particles of type 2. The energy En 
denotes the total (kinetic plus potential) energy of type- 1 
atoms in the i th direction. Note that our current interest 
in xi implies the ensemble average in Eq.@ is taken only 
over the distribution function for type-1 atoms. Based 
on the results of Fig. |2 we can assume equal numbers of 
atoms Ni = N2 without loss of generality. 

Immediately after a collision, only the kinetic energy 
(KE) has changed, so that 



A(KE lx - 
- mi A(vi 



KE lz ) 



Since A^i has no position dependence, we can remove 
the type-2 particle density n 2 from the average, yielding 

< Xl ) = 5 TOl ( " 2 ) ai2 ( Urcl A ( U lz ~ V lz)) ■ 

We define CM and relative velocities in the usual man- 
ner, 



VcM = 



mivi + m 2 t>2 



mi + m 2 
Vie} = Vi-V 2 , 



(5) 



which gives 



J lx 



rn 2 



mi + m.2 
m 2 

mi + m 2 



re\x relzy 



— VcMz V-cc\ z ) 



Since the collision leaves Vcm and |tT re i| unchanged, but 
randomly rotates the direction of v re \, we are left with 



1 mi m 2 

2 (mi + m 2 ) 
m 2 



n 2 ) a 12 



x ( v vci 



( v ielx ~ v ie\z) 

mi + m 2 

+ 2 (VqMx Vrelx ~ VcMz ^rolz) ^> 

Here we need only consider the quantities immediately 
before the collision since there is no preferred direction 
for v Te \ after the collision. Note also that ( Vcmi VreU ) 7^ 
for a gas with a cross-dimensional energy anisotropy. 

Calculating these ensemble averages for arbitrary 
masses and energy anisotropies (which are in general dif- 
ferent between species) is non-trivial, but some simple 
approximations may be used. For the first term it is easy 
to show that for Boltzmann distributions in equilibrium 
one has 



(6) 



For cross-terms of the form ( VcMx ^rels ) 1 which vanish 
under equilibrium conditions, one can still consider the 
limit of small deviation from thermal equilibrium. In this 
case we obtain, similar to the above result, 



— - ( Vrcl ) ( VcMx Vrdx , 



We verified by analytic means and with Monte Carlo in- 
tegrations of Gaussian distributions that these approxi- 
mations are reasonable for small anisotropies. 
Combining these results gives 



mi m 2 



m 2 



--(r 12 

3 mi + m 2 \ mi + m 2 

+ 2 (VcMx Vrclx — VcMz V Ie l z ) 
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The collision rate Ti2 describes the rate per type-1 par- 
ticle of collisions with particles of type 2. If we finally 
substitute back with v\ and v 2 , and use ( vu u 2 i ) = 0, we 
obtain 



Xi 



2 m 2 
q 1 12 7 



m 2 ) 



(2 mi + m 2 ) xi - Xi } (7) 



where we have dropped the angle brackets now for sim- 
plicity. Note that we have used the fact that the mean 
kinetic and potential energies in a given direction are 
equal in the collisionlcss regime. 

For concreteness, we now associate the type-1 par- 
ticles with the fermions and type-2 particles with the 
bosons. The time-dependence of X2, defined in analogy 
with Eq.©, is obtained by swapping 1 <-> 2 in Eq.© 
and adding a term describing the effect of boson-boson 
collisions. This term was calculated from Enskog's equa- 
tion in Ref.0, giving a = 5/2 in the limit of small energy 
anisotropy. This yields the final result, 

d 2 m 2 . . . , 

-j- Xi = -o7 i (2mi+m 2 )xi-rniX2 

dr S [mi + m2) 

d 2 mi . 

~T X2 = -o? ; Y2 i( 2m 2 + mi)x2-m 2 xij 

a t 6 (mi + m 2 j z 

--X2 , (8) 
a 

where we have introduced the dimensionless time 
t = t and the ratio of collision rates 7 = T 22 /ri 2 
(note that Ti 2 = T 2 i for Boltzmann distributions with 
Ni = N2 under conditions of thermal equilibrium). 

Equation © is the main result of our analysis. The 
results of exponential fits to the time-evolution of Eq. © 
are shown as the open points in Fig. [3] The slight over- 
estimation (underestimation) for light (heavy) fermions 
is a result of our approximation in Eq. © . Although the 
full solution to Eq.© requires prior knowledge of the 
inter-species cross-section <7i 2 (through 7), we now show 
that the solutions in the limits 7 = and 7 — > 00 pro- 
vide tight bounds on j3. These two limits are shown as 
the solid lines in Fig. |3J 

In the case 7>1, where boson-boson collisions occur 
much more frequently than boson-fermion collisions, we 
can assume the energy of the bosons is always isotropic 
and take X2 = for all times. In this limit, Eq.© is 
easily solved, giving 



Xi(t) = Xi(0)exp 



with 0i given by 



2 1 



(9) 



The solution to Eq.© in the limit 7^0 does not cor- 
respond exactly to a simple exponential but rather to 
the sum (or difference) of decaying exponentials. Heavy 
fermions (77 > 1/2) rethermalize slowly at first and 



then more quickly as the energy of the bosons becomes 
isotropic. In this case, we obtain an upper limit on 
f3 by simply taking the small-r expansion of the solu- 
tion to Eq.©, giving (3 U = (3/2) (1 - t/T 1 . For light 
fermions (77 < 1/2), we obtain j3 u by fitting the solution 
of Eq.© to the simple exponential decay minimizing the 
chi-squared error. Our values for (3 U and fi\ are shown 
in Fig. 13 as the upper and lower solid lines, respectively. 
The results from the Monte Carlo simulations fall within 
these limits as expected. Simulations with different atom 
numbers and cross-sections gave similar results. 



IV. DISCUSSION 

Before proceeding, we first consider the effect of the 
finite size of the initial perturbation. Our analysis has 
so far been limited to the case of arbitrarily small initial 
anisotropies. As discussed in the introduction, however, 
experiments (as well as Monte Carlo simulations) require 
large anisotropies for good signal-to-noise ratios. It was 
noted in Ref. 5] that the mean collision rate per particle 
( r co n ) = (n) a (v m \) actually changes during the relax- 
ation process, because of the redistribution of energy. In 
Eq.© we have defined the relaxation rate in terms of 
the final equilibrium collision rate. Because of the ex- 
ponential nature of the relaxation, however, the bulk of 
rethermalization occurs at initial times. We therefore ex- 
pect that the observed number of collisions required for 
equilibration will be approximately rescaled by a factor 
of (T co \i(t — > 00) )/( r co n(t = 0) ). As shown in Appendix 
1X1 this effect can be accounted for by multiplying j3 by 



A(fi) = n 



3/2 



tan -1 



-(1 + 20) 
, -1 



(10) 



The function A(O) goes to zero at cither limit of fl 
(0 or 00) and reaches a smooth maximum of 1 at Q = 1. 

To test our prediction for A(f2) given by Ea. (jl()|l . we 
performed a series of Monte Carlo simulations, both for 
single-species relaxation of 87 Rb and for relaxation of 40 K 
atoms in the presence of 87 Rb, as shown in Fig. 0J The 
model and the Monte Carlo show excellent agreement 
over a wide range of Q. 

Finally, we predict in Table ^ values of j3 for various 
mixtures of bosonic and fermionic alkali atoms. The 
nominal value is the average of /3 U and (3i and the un- 
certainty is half the difference. To date, all Bose-Fermi 
mixtures produced in experiments have light fermions 
(7/ < 1/2) , where the maximum uncertainty in our predic- 
tion is about 20%. Since measurements of the magnitude 
of the scattering length depend on J|, the model 

will introduce only a small uncertainty in any determi- 
nation of IobfI as described in this work. As a test of 
the performance of the prediction, we consider the 40 K- 
87 Rb system, for which our model predicts (3 = 2.1 ± 0.4. 
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FIG. 4: Ratio of rethermalization and collision rates as a func- 
tion of the initial perturbation, (a) Single-species relaxation 
of 87 Rb. The points are results from the Monte Carlo simula- 
tions, and the solid line is the best fit to Ea. llOH . which gives 
a = 2.64 ± 0.07 for the limit of small initial anisotropy. (b) 
Number of collisions per fermion ( 40 K) for rethermalization 
with 87 Rb. Points are from the Monte Carlo simulations, and 
the upper and lower dashed lines combine Eq. llOU with p u 
and Pi, respectively. The solid line is the best fit to Eq. dlOL 
giving f3 = 2.26 ± 0.09 for vanishing initial anisotropy. 



Monte Carlo simulations over a wide range of Nb,Nf, 
and n gave (3 = 2.23 ± 0.07 (in the limit ft = 1), in 
excellent agreement with the prediction. 



V. CONCLUSION 

We have investigated the cross-dimensional rethermal- 
ization of fcrmions in the presence of bosons, focusing 
attention on the effects of the mass difference between 
species and the finite departure from equilibrium nec- 
essary for experiments. Monte Carlo simulations were 
performed over a wide range of parameters, and a simple 





6 Li 


40 K 


7 Li 


2.4 ±0.5 


8±2 


23 Na 


f.9±0.3 


3.3 ±0.8 




1.8 ±0.2 


2.5 ±0.5 


87 Rb 


1.64 ±0.14 


2.1 ±0.4 


133 Cs 


1.60 ±0.09 


1.9 ±0.4 



TABLE I: Predicted values of /3 for Bose-Fermi mixtures of 
commonly used alkali atoms. The nominal value is the aver- 
age of p u and Pi, and the uncertainty is equal to half of the 
difference p u — Pi- 



analysis based on Enskog's equation was developed that 
reproduces the results from the simulations. We have fur- 
ther used the model to predict the number of collisions 
per fermion needed for rethermalization for a variety of 
alkali Bose-Fermi mixtures. 

In this work, we have restricted ourselves to rethermal- 
ization due to energy-independent s-wave collisions, but 
the Monte Carlo simulations can easily be extended to 
include p-wave collisions Q , resonant scattering 0] , or 
damping of more complicated collective excitations 0. 
In addition, the simulations have allowed us to investi- 
gate various effects not accounted for in our model, such 
as the differential gravitational sag for species with dif- 
ferent masses 01 , an d the possibility of slightly different 
initial energies and anisotropies between species. The 
simulations could be trivially extended to Fermi-Fermi 
or Bose-Bose mixtures 0] and could accommodate 
in a straightforward way relaxation in the hydrodynamic 
regime, Bose- Einstein and Fermi-Dirac distributions [l5| . 
heating and loss during the rethermalization, and the ef- 
fects of anharmonic trapping potentials. 



APPENDIX A: CALCULATION OF A(O) 

As discussed in Sec. IIVI the rescaling A(S1) is given by 
the ratio of the final and initial collision rates, 



A(O) 



(n B (t — > OO)) (Vreljt — > OO)) 
(n B (t = 0)j («rel(* = 0)) 



We assume here that the mean energy in a given direction 
is the same for fermions and bosons at the beginning and 
end of the relaxation, and denote these energies 



(E x ,E y ,E z 
kB ^ 



1 ±2fi 



(n,n,i) , t = o 



(Al) 



(1,1,1) 



t 



OO 



A simple calculation for Gaussian distributions gives 
the mean boson density (averaged over the fermion dis- 
tribution function), 



(n B ) cx (E x E y E z 



-1/2 
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which, using Eq. (|Al|l . gives 

(n B (oo)) „ / 3 



.3/2 



(n B (0)) ""Vl + 20 
For (i» re i) one finds for our initial conditions, 

fl tan -1 - 1 



Note that for £2 < 1 the fraction in the square brackets 
can be written (1 — fi) -1 / 2 tanh -1 (1 — f?) 1 / 2 , which is 
again real- valued. The product of Eas. lA2(l and (|A3(I 

(A2) 

was presented as A(f2) in Eq. l|l(J|l . 



(«rcl(0)) 



' 3 2 fc B r 

1 + 20 7T/Z 



1 



where fi = mp ms / (mp + tub) is the reduced mass, as 
in the text. Comparing to the equilibrium value gives 



(^rcl(oo)) 
(Urel(O)) 



-(1 + 212) 



1 + Sl 



tan- 1 VTi~l 
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